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FIELD OF THE INVENTION 

The present invention relates to a method of generating a hybrid grid allowing 
modelling of a heterogeneous formation crossed by one or more pipes. 

The method is more particularly applied to formation of a grid suited to an 
5 underground reservoir crossed by one or more wells, in order to model displacements of 
fluids such as hydrocarbons. 

BACKGROUND OF THE INVENTION 

Grid generation is a crucial element for the new generation of reservoir simulators. 
Grids allow to describe the geometry of the geologic structure studied by means of a 

10 representation in discrete elements wherein simulation is performed according to a 
suitable numerical pattern. Better comprehension of physical phenomena requires 3D 
simulation of the multiphase flows in increasingly complex geologic structures, in the 
vicinity of several types of singularities such as stratifications, faults, pinchouts, 
channels and complex wells. All this complexity has to be taken into account first by 

15 the grid which has to reproduce as accurately as possible the geologic information in its 
heterogeneous nature. 

Grid modelling has made great advances during the past few years in other fields 
such as aeronautics, combustion in engines, structure mechanics, etc. However, the 
gridding techniques used in the other fields cannot be applied as they are to the 
20 petroleum sphere because the professional constraints are not the same. For example, in 
reservoir simulation, the numerical patterns are constructed from control volumes in 
order to better respect the mass conservation in the case of transport equations of 
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hyperbolic nature. The grid must be a « block-centered » type grid, i.e. the nodes must 
be inside each layer and the boundaries of each block must follow the interface between 
the layers. Now, if this constraint was not taken into account, the nodes would naturally 
be placed along the faults and along the stratification boundaries. The consequence of 
5 this would be that these interfaces would pass through the control volume used. The 
saturation, constant in the control volume, could not respect then the discontinuity and 
the results would not be accurate. It is therefore necessary to develop new techniques 
that are better suited to the petroleum sphere requirements. 

Cartesian grids, which are commonly used in current commercial simulators, are 
10 unsuited for solving these new problems posed by the development of petroleum 
reservoirs. Cartesian grids, based on parallelepipedic elements, do not allow 
representation of such complex geometries. 

There is a well-known method of generating structured 3D hexahedral grids of CPG 
(Corner-Point-Geometry) type which respects the geometry of the bodies. It is described 
15 in patent FR-2,747,490 (US-5,844,564) filed by the applicant and also in the following 
publication : 

- Bennis Ch. Et al. « One More Step in Gocad Stratigraphic Grid Generation » : 
Taking into Account Faults and Pinchouts ; SPE 35526, Stavanger, 1996 . 

This grid type is more flexible than cartesian grids because it consists of any 
20 hexahedral elements that can be degenerated. It strictly respects the horizons, the faults 
and it allows to represent certain unconformities such as pinchouts because its 
construction is based on these elements. However, this type of grid does not allow to 
solve all the geometric complexities such as, for example, circular radial grids around 



complex wells. It is possible to form separately the grid of the reservoir and the grids 
around the wells but it is difficult to represent several objects in the same CPG type 
reservoir grid because of connection problems linked with the structured nature of the 
grid. 

Another approach is also known where 3D grids only based on tetrahedral Delaunay 
elements, with a circular radial refinement around the wells, are automatically 
generated. The advantage of such an approach is that it is entirely automatic and does 
practically not require the user's attention. However, this method has drawbacks which 
make the results obtained difficult to use : 

- there are on average five times as many grid cells as in a CPG type grid for the same 
structure, which is very disadvantageous for simulation calculations, 

- unlike the structured grids which are easy to visualize, to explore from the inside 
and to locally modify interactively, it is very difficult and sometimes impossible to 
properly control the tetrahedral grids because of their size and especially because of 
their non-structured nature. This poses problems for validating the grid from a 
geometric point of view as well as for understanding and validating the result of a 
simulation on this type of grid. 

Other approaches are also well-known, which allow to generate grids, notably grids 
based on control volumes generated from a triangulation (Voronoi and CVFE), 
associated with techniques of aggregation of the triangles (or tetrahedrons) into 
quadrangles allowing the number of grid cells to be reduced. Although promising 
results were obtained with these new grids, precise representation of the geologic 
complexity of reservoirs and wells remains a subject for research and development. In 
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fact, these grids are rather 2.5D (i.e. vertically projected) and their 3D extension appears 
to be very difficult. Despite their hybrid aspect, they remain entirely unstructured and 
would therefore be very difficult to manage and to handle in real 3D. Furthermore, 
taking account of real 3D faults and deviated wells would greatly increase this 
5 difficulty. 

SUMMARY OF THE INVENTION 

The method according to the invention allows to generate a hybrid grid suited to a 
heterogeneous formation crossed by at least one pipe of known geometry (such as an 
underground reservoir crossed by one or more wells), in order to form a model 
10 representative of fluid flows in this formation in accordance with a defined numerical 
pattern, the structure of the formation being known a priori from available data acquired 
by means of in-situ measurements, analyses and/or interpretations of formation images 
(seismic images for example, in the case of a reservoir). 

The method is characterized in that it comprises associating a first structured grid 
15 for gridding of the formation while respecting the discontinuities thereof with second 
structured, radial type grids for gridding of the zones around the wells, these second 
grids allowing to respect constraints linked with the flows in the wells, and non- 
structured transition grids between the first grid associated with the formation and the 
second grids associated with the wells. 

20 Gridding of the heterogeneous medium is obtained for example by importing each 

second structured grid into a cavity formed in the first structured grid, the size of this 
cavity being sufficient to allow formation of a non-structured transition grid between 
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the first structured grid associated with the formation and the second structured grid 
associated with each well. 

The non-structured transition grid based on any polyhedrons or canonical 
polyhedrons such as pentahedrons, tetrahedrons, pyramids, etc., can be formed by 
5 respecting constraints linked with the numerical pattern. 

The non-structured transition grids are advantageously modelled with the structured 
well grids by applying a formalism known in the art, referred to as « generalized map » 
formalism, the grid of the formation being structured matrically, globally or in faulted 
blocks. 

10 The global hybrid grid is thus obtained by combination of several grid types : a 

structured reservoir grid, a radial grid around each well, also structured, and non- 
structured transition grids which connect the previous two grid types. Each one of these 
grids has its own formation, representation and exploration methods. The structured 
aspect is thus degraded only at the points where this is strictly necessary. This « object » 

15 approach affords both the advantage of structured grids for control and comprehension 
of the reservoir and the flexibility of non-structured grids in complex zones. Complexity 
is introduced only where it is strictly necessary. The independence of these gridding 
modes therefore allows separate extraction, management and representation of the well 
grids and of the interstitial grids included in the reservoir grid. 

20 Using a reservoir simulator of a well-known type, such as ATHOS™ or SCORE™ 

for example, for a reservoir provided with a hybrid grid obtained by means of the 
method, allows production simulations to be performed. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Other features and advantages of the method according to the invention will be 
clear from reading the description hereafter of non limitative examples, with reference 
to the accompanying drawings wherein : 

5 - Figure 1 shows a diagrammatic example of a hybrid grid of a reservoir crossed by 
two wells, consisting of a first structured grid for the reservoir, a second structured grid 
for the zones around the wells and transition grids between the two grid types, 
Figure 2 shows an example of a structured grid of a faulted reservoir, 
Figure 3 shows an example of a radial grid around a vertical well, 
10 - Figure 4 shows an example of grid of a horizontal well, 

- Figure 5 shows, in 2.5D, an example of a gridded reservoir where cavities are 
provided for gridded wells, before the stage of creation of non-structured interstitial 
grids intended to connect them together, 

- Figure 6 shows five wells provided each with a radial grid, included in a gridded 
15 reservoir, by means of non- structured transition grids based on any polyhedral grid 

cells, 

- Figure 7 shows separately the structured reservoir grid with the cavities provided 
therein in order to include the additional elements : gridded wells and interstitial grids, 

- Figures 8A, 8B, 8C, 8D show elementary well grids visualized individually 
20 according to different modes, 
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Figures 9A, 9B, 9C, 9D show various elementary transition grids, visualized 
individually according to different modes, also, allowing their integration into the 
reservoir grid, 

- Figure 10 shows an example of a hybrid reservoir grid with transition grids 
5 consisting for example of pentahedrons, between the reservoir grid and several well 

grids, 

Figure 1 1 shows a model for the matrical representation of a structured grid around 
a well, 

- Figure 12 is a graphic representation of a connection between strands, within the 
10 scope of the modelling technique referred to as « generalized maps » technique used for 

generating non -structured grids, 

- Figures 13 A, 13B are graphic representations of connections by means of simple 
arcs, and 

- Figures 14A, 14B are graphic representations of connections by means of double or 
15 triple arcs respectively. 

DETAILED DESCRIPTION 

Modelling of the reservoir is obtained by combining elementary grids of different 
types. Each elementary grid is considered to be a full object with its own data model, its 
own generation methods and its own representation methods. Generation is carried out 
20 in stages with addition/subtraction of grids. 

1) In order to represent the reservoir as a whole, an i, j, k structured grid of a type 
known to specialists, referred to as CPG (Corner Point Geometry), as described in the 



aforementioned patent FR-2,747,490 is used for example. The reservoir can be faulted 
with downcreep of a block in relation to the other. The major horizons and faults are 
first modelled by continuous surfaces from data resulting from an interpretation of 
seismic images of the subsoil or from data obtained during drilling (well markers). The 

5 geologic structure is then divided into faulted blocks resting on these surfaces. These 
blocks are individually gridded, then reassembled. Gridding of a block first consists in 
gridding the edge surfaces, then the inside is populated by transfinite interpolation of 
the edge surface grids. Relaxation techniques are then applied to the edge surfaces and 
to the inside so as to harmonize and to regulate the grid. The grid thus obtained strictly 

10 respects the horizons, the faults and it allows to represent certain unconformities such as 
pinchouts. It meets all the constraints of geologic nature. 

2) A well trajectory is drawn synthetically or imported. A structured radial grid is 
then generated around each well in order to take account of the particular constraints 
linked with the flows in the vicinity of these wells. 

15 In the example shown in Fig.3, the structure grid around a vertical well is of circular 

radial type. It is also a CPG type grid. Its generation first consists in sampling a disc at 
r, 0 in the horizontal plane. The 2D grid thus obtained is then projected vertically upon 
the various layers of the reservoir grid. Here, the i, j, k of the matrical structure 
correspond to the samplings at r 5 9 and z respectively (see Fig.l 1). 

20 The grid around a horizontal well (Fig. 4) is i, j, k structured, it is of the same type as 

that of the reservoir, except that a well cannot be faulted. It is also obtained by 
projecting vertically upon the various layers of the reservoir grid a 2D grid belonging to 
a horizontal plane. 
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3) This radial grid is then inserted around the or around each well in the global 
reservoir grid. A cavity is therefore first created in the reservoir grid by deactivating all 
the grid cells in contact with well grid cells (Figs.5, 6). The space freed between the 
reservoir grid and the well grid must be sufficient to allow convenient formation of a 

5 transition grid. It can represent for example about the equivalent of two grid cell layers. 

4) A non- structured transition grid is then generated in this cavity (Fig.7) in order 
to connect the structured radial grid around the well to that of the reservoir best 
respecting the constraints linked with the numerical pattern. The user can deactivate the 
grid of a well any time by reactivating the grid cells of the corresponding cavity in the 

10 reservoir grid. 

The transition grid can for example consist of polyhedrons with any number of 
sides or canonical polyhedrons (tetrahedrons, pentahedrons, pyramids, etc.) according to 
the numerical pattern used, without the overall hybrid approach proposed being 
affected. 

15 Example of modelling of a hybrid grid 

The reservoir grid and each well grid are modelled, for each faulted block of the 
reservoir, by matrical structures of points or cells comprising each eight points. Because 
of the structured nature of the grids, the topological links between the various grid cells 
are implicitly contained in the matrical structure. 

20 Transition grids are more difficult to manage because of their non-structured nature 

and because they can contain polyhedral grid cells whose number of sides varies from 
one cell to the other. 



An advantageous solution for facilitating management of this new grid type, 
allowing to browse it and to surf it efficiently, consists in using the topological model 
referred to as « generalized maps or G-maps ». This model known to specialists is for 
example described by : 

5 - Edmond J. : « A Combinatorial Representation for Polyhedral Surfaces », Notice 
Amer. Math. Soc, 7, 1960, or by : 

- Fortune S., 1992 : Voronoi diagrams and Delaunay triangulations, pp.225-265 of 
D.Z. Du & F.K. Hwang (eds.), Computing in Euclidean Geometry, 2 nd edn. Lecture 
Notes Series on Computing, vol.4, Singapore, World Scientific. 

10 Generalized maps are based on a formal algebraic approach that is briefly reminded 

hereafter. 

In 3D, the elements which constitute a generalized map are (D, cto, a u 0C2, CX3), 
where D is a finite set of elements called strands, and elements {0^} are involution on D 
type functions, associating the strands two by two at most, which are therefore 
15 conveniently referred to as links. Figs. 11 to 14 show concrete geometric representation 
examples. Link Oo is in the form of a dotted segment (Fig. 11) and links oti, 0C2 and 0C3 in 
the form of arcs, respectively simple (Fig. 12), double (Fig. 13) and triple (Fig. 14). 

According to another known approach, generalized maps are considered as graphs 
whose strands form the nodes and the links form the arcs : link Oo between two strands 
20 can be used for representing the edge of a side, links a, for connecting two edges of a 
side, links 02 for linking two sides of a cell together and links oc 3 for sticking two cells 
together. 
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This model of generalized maps involves a small number of formal objects and an 
operation which, by associating additional information with a topology, allows to locate 
the objects defined in space and to account for their appearance, which is referred to as 
plunge and, in the present case, plunge in a 3D space. 

5 It affords the advantage of being independent of the dimension of the objects. All 

the objects can be represented with the same data structure and handled with the same 
methods. This approach makes it possible to handle objects created with heterogeneous 
topological models. It is therefore well-suited for implementing the method according to 
the invention with its stage of creation of a non-structured grid linking two different 

10 structured grids together. 

The generalized map concept for modelling the transition grid is applied by creating 
a certain number of objects of different types which refer to one another. These objects 
materialize the topological network and its various plunges in a 3D space. Concretely, 
in order to allow browsing the grid, a topological network is constructed parallel to the 
15 geometric data commonly handled in a grid, the points, the edges, the sides and the 
cells. Furthermore, crossed links are established between the topological network and 
geometric data. 

Objects 

The various objects handled within the scope of the application performed here of 
20 the generalized map method are as follows : 

1) The Transition Grid object which contains all the topology, the geometry and 
the physical data. It consists of a GMap type object which represents the topological 
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network and of a Plunge type object which materializes the plunge in the physical world 
according to our application. 

2) The GMap object 

The topological model is entirely contained in a graph consisting of a list of Strands 
5 connected to one another. Any operation performed on the generalized map amounts to 
an operation on the Strands network. The GMap object type has methods allowing easy 
circulation in the topological network representing the grid, i.e. to go from one Strand to 
another. 

3) The Strand object 

10 Each Strand is defined by four references to other Strands (corresponding to links 

Oo, oti and a 3 ) and by four other references to the plunge in the 3D space, notably at a 
Point, an Edge, a Side and at a Cell to which it is connected. 

4) The Plunge object 

It is defined by four lists : 
15 - a list of Points (the grid points), it is the plunge of dimension 0 of the GMap, 

- a list of Edges (plunge of dimension 1 of the GMap), 

- a list of Sides (plunge of dimension 2 of the GMap), and 

- a list of Cells (plunge of dimension 3 of the GMap). 

The Plunge object also contains its own methods of creating and handling the data it 
20 contains according to the use that is made thereof. Besides, the GMap is created from its 
Plunge. 
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5) The Point object 

A Point is defined by its coordinates x, y, z and by a list of attributes, notably scalar 
or real petrophysical values that are associated therewith. 

6) The Edge object 

5 It is defined by a reference in the GMap to a Strand which represents an end of the 

Edge. This gives a preferential access to the topological network and simultaneously 
allows to go from the plunge to the strands graph. For example, link Oo of the Strand in 
question leads to the Strand representing the other end of the Edge. It is furthermore 
defined by a list of attributes, notably scalar or real petrophysical values that are 

10 associated therewith. 

7) The Side object 

This type of object allows to handle directly the interfaces between the cells as well 
as the outer sides of the grid. A side is defined by a reference in the GMap to a Strand 
which represents a vertex of the edge polygon of the Side. This also gives a preferential 
15 access to the topological network. The Strands representing the other vertices of the 
polygon are accessible by iterative applications of the relation o<,oai by starting from 
the initial Strand and eventually coming back to this same Strand. It is also defined by a 
list of attributes relative to the Side (for example, scalar or vector petrophysical values). 

8) The Cell object 

20 The type of Cell object is defined by a list of references in the GMap to Strands, 

each one representing a half-Side of the Cell. This allows access to the topological 
network from the Cells. It is also defined by a site (coordinates of the center of mass of 
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the Cell) and by a list of attributes specific to the Cell or to its site (scalar or vector 
petrophysical values for example). 

Graphical representation and exploration 

A graphical representation is a very efficient and even essential means for 
5 controlling and validating the construction of a grid and the simulation results. 
Concerning the construction, the geometry of the grid generated is generally first 
visually controlled. If this is not sufficient, local or global quality criteria with which 
statistics are established can be calculated and visualized on the grid by means of a 
colour scale. Flow simulation consists in calculating the variations with time of certain 
10 petrophysical parameters by taking account of the hypotheses that initially condition the 
flows. Simulation validation also involves visualization of these parameters on the grid 
(preferably by means of a colour scale). As the grids concerned are 3D grids, tools 
allowing to explore the grid from the inside by visual browsing are required. Graphical 
representation and browsing in the grids, presented hereafter, are a good illustration of 
15 the flexibility and the modularity of the hybrid approach proposed and of the efficiency 
and the adequacy of the data model selected. 

The hybrid grid, considered as a set of independent entities : the elementary grids, is 
constantly visualized in a main window. The user can select at any time an elementary 
grid and visualize it with its specific methods in a secondary window which contains 
20 only the elementary grid selected. Actions on the elementary grid have automatically 
repercussions on the entire hybrid grid visualized in the main window. An elementary 
grid can thus be visualized and explored as a full entity and it can be viewed in the 
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global context. The visualization methods differ according to whether the elementary 
grid is structured (reservoir grid and well grids) or not (transition grids). 

Examples of functionalities specific to structured grids 

In the case of a structured grid (reservoir and well), visualization is simple and 
conventional. It consists in two main functionalities : 

visualization of the external envelope of the grid with the possibility of peeling it in 
the 3 directions i, j, k separately, 

- simultaneous or separate visualization of three matrical cell slices i=cste, j=cste and 
k=cste, with the possibility of moving them in the block. 

Examples of functionalities specific to non-structured grids 

In the case of non-structured transition grids, other, more elaborate visualization 
modes are preferably used. Five functionalities are mainly used : 

visualization of the external envelope with the possibility of concentric peeling, 
topologically speaking, 

- visualization of the cells crossed by a cutting plane orthogonal to an axis of 
coordinates x, y, z or any axis, 

- visualization of the trace of the cells on the cutting plane, 

- visualization of the grid sites when they are intially given and 
visualization of the grid cells in full or scattered mode. 

Of course, for the two grid types, it is possible to visualize a property or scalar value 
by means of a colour scale. 
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All these functionalities require easy and optimum viewing of the non-structured 
grid. This is possible by using the formalism referred to as generalized map formalism. 

Figs. 7 to 10 clearly illustrate the potential afforded by the hybrid grid method 
proposed, i.e. harmonious integration of a structured grid following a topological model 
(well grid) into another structured grid (reservoir grid) following a different topological 
model, by means of a non- structured transition grid. The independence of these models 
therefore allows extraction and separate representation of the well grids and of the 
interstitial grids included in the reservoir grid in order to represent, handle and explore 
this type of data. 



